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CRACK  TIP  FIELDS  FOR  FAST  FRACTURE 
OF  AM  ELASTIC-PLASTIC  MATERIAL 


by 

J.  D.  Achenbach,  M.  F.  Kannlnen,  and  C.  B.  Popelar 


INTRODUCTION 


An  expanding  interest  la  developing  in  the  analysia  of  rapid 
unstable  crack  propagation  in  which  dynamic  (inertia)  effects  may  not  be 
negligible.  While  this  Interest  has  been  motivated  in  part  by  intellectual 
curiosity  in  a  new  and  fertile  applied  mechanics  research  area,  there  is 
also  a  growing  realization  that  many  engineering  applications  for  this 
technology  exist.  There  are  many  practical  situations  in  which  the  result 
of  large  scale  unstable  crack  growth  is  completely  unacceptable.  To  pre¬ 
clude  such  catastropic  occurrences,  it  is  Important  to  develop  treatments 
for  the  arrest  of  rapid  crack  propagation.  Because  crack  arrest  must 
logically  be  viewed  as  the  termination  of  fast  fracture,  attention  must 
be  focused  on  the  dynamic  crack  propagation  process. 

For  the  most  part,  dynamic  fracture  mechanics  is  now  based  upon 
elastodynamic  solutions;  for  example,  see  Kannlnen  (1978).  Yet,  most  struc¬ 
tural  components  where  fracture  is  a  concern  employ  tough  ductile  materials 
where  it  is  unlikely  that  the  basic  assumptions  of  linear  elastic  fracture 
mechanics  (LEFM)  are  valid.  In  addition,  recent  work  has  begun  to  raise 
serious  questions  about  the  extent  to  which  the  dynamic  fracture  toughness 
is  truly  a  material  property,  even  when  LEFM  conditions  are  satisfied. 

Thus,  it  is  necessary  to  develop  nonlinear  dynamic  elastic-plastic  frac¬ 
ture  mechanics  treatments. 


2 


Because  unloading  may  take  place  In  dynamic  crack  propagation, 
a  da formation  theory  formulation  is  undesirable.  A  flow  theory 
approach  is  required  to  account  for  the  elastic  unloading  in  the  wake  of 
the  crack  tip.  An  asymptotic  solution  near  the  tip  of  a  crack  propagating 
dynamically  in  antiplane  strain  (Mode  III)  conditions  using  such  an  approach 
has  already  been  given  by  Achenbach  and  Kaxmlnen  (1978).  Their  results 
determined  the  form  of  the  ersek  tip  singularity  and  the  angular  position 
at  which  elastic  unloading  conmences.  Clearly,  this  kind  of  information 
is  needed  in  order  to  devise  a  finite  element  or  other  solution  procedure 
for  dynamic  crack  propagation  and  arrest  analyses. 

The  work  described  in  this  paper  extends  the  Achenbach-Kanninen 
approach  to  Mode  I  crack  growth.  Results  for  both  plane  stress  and  plane 
strain  conditions  are  given.  These  solutions  are  based  upon  J2  flow  theory 
with  a  bilinear  stress-strain  relation  which  allows  elastic  unloading. 

The  problem  formulation  leads  to  a  nonlinear  eigenvalue  problem  which  is 
solved  using  an  iterative  numerical  solution  procedure.  Confidence  in 
the  results  was  attained  by  comparisons  with  the  results  obtained  by 
Amazlgo  and  Hutchison  (1977)  for  the  special  case  of  quasi-static  crack 
growth  at  low  crack-tip  speeds  when  the  dynamic  effects  are  negligible. 
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THE  SOLUTION  PROCEDURE 


Problem  Formulation 

The  fields  of  stress  end  deformation  ere  referred  to  e  coordinate 
system  whose  origin  is  attached  to  the  moving  crack  tip.  The  system  of  coor¬ 
dinates  is  shown  in  Figure  1.  The  crack  is  located  in  the  (x^z^)-plane  where 
the  axis  coincides  with  the  crack  front  and  x^  is  the  direction  of  crack  * 
advance.  The  relevant  displacement  components  are  (x^.x^.t)  and  ^(x^^tt), 
where  t  is  time. 

In  the  following  the  material  derivatives  with  respect  to  time  will 
frequently  be  needed.  These  are  defined  as 


is : 


+  [v(t)l 


2  jr 


(1) 

(2) 


where  v(t)  is  the  speed  of  the  crack  tip.  Notice  that  v  need  not  be  constant 
and  is  only  subject  to  the  conditions  that  v(t)  and  dv/dt  are  continuous  functions. 

For  the  case  of  plane  stress,  the  non-vanishing  stress  components  are 
Ojj,  o12  »  o21,  and  The  equations  of  motion  are  of  the  form 

o  .  .  -  pu  *  Y»5  “  1*2  (3) 

ye, 6  y 

where  the  second-order  material  time-derivative  is  defined  by  (2)*. 

The  constitutive  equations  in  the  elastic-plastic  material  take  into 

account  strain-hardening  characterized  by  Jj  flow  theory  and  a  bilinear  effective 

*ln  the  following,  Greek  minuscules  have  the  values  1,2  whereas  Latin  minuscules 
have  the  values  1,2,3. 
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stress-strain  curve.  This  curve  is  shown  in  Figure  2.  Let  oq  denote  the 

yield  stress  in  tension,  E  Young's  modulus,  end  a  “  E^/E,  where  Et  is  the 

slope  of  Che  bilinear  stress-strain  relation  in  tension  for  stresses  in  excess 

of  o  .  The  effective  stress  is  defined  as 
o 


a 

e 


*ij  *ij^ 


1/2 


(4) 


where  s^  Is  the  stress  deviator.  The  constitutive  relations  for  an  elasti¬ 
cally  isotropic  solid  may  then  be  written  following  Amazlgo  and  Hutchinson 
(1977) : 

loading  (o#  >  o): 

Et  *ij  "  a*(1  +  v)6ij  “  ^kk^ij1  +  <|  °e>(1  -  (5) 

unloading  (a  <  o): 

© 

Et®ij  "  «[(1  +  v)6±j  -  (6) 

where  v  is  Poisson's  ratio. 

Field  Equations  for  a  Near-Tip  Analysis 

In  this  paper,  certain  a-priori  assumptions  are  made  with  respect  to 
the  general  nature  of  the  deformation  in  the  immediate  vicinity  of  the  crack  tip. 
Explicit  expressions  for  the  near-tip  fields  will  be  derived  using  these  assump¬ 
tions.  It  will  then  be  shown  that  the  required  boundary  and  continuity  condi¬ 
tions  are  satisfied. 

It  is  assumed  that  ahead  of  the  crack  tip  a  plastic  loading  zone 
exists  which  is  bounded  in  the  near-tip  region  by  radial  lines  emanating  from 
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the  crack  tip  at  angles  8  ■  +0p.  This  Is  shown  in  Figure  1.  An  elastic  unload¬ 
ing  rone  Is  assisted  for  |  &|  >  8^.  The  possible  presence  of  a  plastic  reloading 
sone  for  | 8|  near  *  is  neglected  in  this  paper.  Conceptually,  it  is  possible  to 
Include  such  a  reloading  sons.  But,  just  as  for  the  quasi-static  case  investi¬ 
gated  by  Amazigo  and  Hutchinson  (1977),  it  is  expected  that  the  Influence  of  re¬ 
loading  is  not  significant  enough  to  justify  the  computational  complications 
that  would  be  required  to  treat  it. 

In  the  plastic  loading  rone,  asymptotic  solutions  near  the  crack  tip 
of  the  general  form 

uy  -  HvBy(8)  r8  (7) 


are  sought.  Here,  K  is  an  amplitude  factor,  while  6^(8)  and  s  are  to  be  deter¬ 
mined.  In  an  asymptotic  analysis,  only  the  lowest  orders  in  r  need  to  be  re¬ 
tained.  This  means  that  3/3t  can  be  neglected  compared  to  -vCOS/Sx^^  in  (1). 
Thus 

(*)  *  -  v(t)  ■£-  (8) 


By  using  the  relations 


and 


CCS  u  - 


sln8  3 
r  38 


sin8I?  + 


cos  8  3 
r  38 


(9) 


(10) 


the  strain  rates  corresponding  to  (7)  are  then  computed  as 
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-  K  v  [s  cose  -  sine]  rB  1 
t22  -  K  v  [s  02  cose  +  *•  cose]  r8"1 

fc12  •  K  v  [s  sine  +  cose  +  a  02  cose  -  V'2  sin6]  r8-1 
where  (  )'  »  d/d 8.  Also  by  definition 

{olj’  °e*  8ij}  "  K  E  {rij(e)»  Ze(e)*  Slj(6)}  r* 

and 

i6±y  6e*  *±jl  -  K  2  V  (1^(8),  Eje),  ^(6)}  r8'1 

By  virtue  of  (8)  and  (9),  the  following  relation  holds 

-  sE^  cose  +  E^j  sine 


Analogous  relations  hold  for  1^(6)  and  (6) .  Also 


Sij  "  lli  ~  3  Skk  6ij 


r  a  (It  R 

e  (2  S1J  Sij, 


By  substituting  (7)  and  the  first  of  (14)  Into  the  equations  of 
(3),  we  obtain  in  the  plastic  loading  zone  near  the  crack-tip 

•  co»  ®  •  sin  6+8  ^21  8in  6  +  J21  C0S  e  "  (“«  COB  8  +  6^  sinC) 

•  cose  -  e*  sme+  s  e22  sine  +  e'2  cose  -  e2  (-s  u,  cose  +  ul  sine) 


where 


B  -  v/c  ,  c  -  (E/p) 


1/2 


(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

motion 

(19) 

(20) 

(21) 
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Substitution  of  (11)  through  (18)  Into  the  constitutive  relations  (5)  and  (6) 
yields 

o[s  cos6  -  sine]  -  o  [(1+v)  iu  -  viy J  +  |  (1-a)  Z"1  Sn  te  (22) 

•  [»  cose]  -  a  [(1+v)  Z22  -  vZjJ  +  (1-a)  Z"1  S22  Z#  (23) 

a|s  sln6  +  .cos6  +  s  U2  cos  6  -  U2  sine]  •  2  a  (1+v)  Z^2 

J  (24) 


+  3  (1-a)  Z"1  Z12  Ze 


The  unknowns  In  these  equations  are  U  and  Z  , . 

y  yS 

The  above  formulation  was  derived  for  plane  stress  conditions.  For 
the  case  of  plane  strain,  we  have  the  condition  ■  0  and  the  nonvanishing 
stress  o^y  It  follows  from  (5)  and  (14)  that  the  additional  relation  required 
for  determining  Z^3  is  # 

“  [(1  +  v)  *33  -  »Hk]  +  T  «  -  »>  h1  S33  K  •  0  <25> 

Ho  other  changes  are  needed. 

Corresponding  equations  for  the  elastic  unloading  region  can  be 
obtained  from  the  work  of  Achenbach  and  Bazant  (1975).  In  the  elastic  region 
the  displacement  rate  must  satisfy  for  plane  strain  the  equation 

<x  +  “)4t.tS  +  ui{,rr  ’  pus"  <26) 

where  X  and  y  are  Lame's  elastic  constants.  We  introduce  the  displacement-rate 
potentials  i  and  f  through 
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i  .  J4.+  2L  .  4  ._y _ hi. 

U1  3*x  3x2  *  U2  3x2  3xx 


It  la  than  not  difficult  to  show  that  and  u2  will  satisfy  (26)  if  4  and  if* 
ara  solutions  of  the  wave  aquations 

♦»rr  "  "7  (*)“  ;  CL  "  (x  +  2w)/p  (28} 

CL 


**ry  "  “7  *  * 

CT 


-  w/p 


In  the  near-tip  region,  the  stresses  In  the  elastic  unloading  zone  nay  be  ex¬ 
pressed  In  tarns  of  4  and  ♦  by 


-vy~1<J22  '  eT  ^  +  ^ 


-vw'^i  *  2  +  (82  -  2) 


where  we  have  used  (8). 


Let  us  now  consider  solutions  of  (28)  and  (29)  of  the  general  forms 


♦  -  k  v  *(eL,e)rp,  eL  -  v(t)/cL 


*  -  K  v  f  (8T,6)rp,  Bt  •  vU)/^ 


Solutions  of  this  kind  follow  Immediately  as  a  slight  generalization  of  Equations 
(22)  and  (23)  of  Achenbach  and  Bazant  (1975)  as 

♦  ($*,6)  -  (1  -  82sin2e)p/2  [A  slnp  (e  -  ir)  +  B  cosp(e  -  *)]  (34) 


9 


¥<e*,0)  -  (1  -  82sin20)p/2  [C  sinp(a>  -  it)  +  D  cosp(u)  -  ir)]  (35) 


where 


Defining 


tan  u  -  (1  -  82)1/2  tan0 

(36) 

tan  e  -  (1  -  62)1/2  tan0 

(37) 

u  «  K  v  U*1^)^1 

(38) 

ve  find  Che  following  relations  on  the  basis  of  (27)  and  (9-10) 
iel 


U*  •  p  cos0  *  -  sine  4’  +  p  sin0  T  +  cosO  ♦* 
6*^  “  P  sin0  4  +  CO8  0  4'  -  p  cos0  ¥  +  sin0T  * 


(39) 

(40) 


The  governing  equations  for  conditions  of  plane  stress  In  the  elastic  region 
can  be  obtained  from  the  above  by  replacing  X  by  v/ (1  -  v) . 

Boundary  Conditions 


The  governing  equations  in  the  plastic  loading  zone,  (19)  through  (24), 
and  the  general  forms  of  the  solutions  in  the  elastic  unloading  zone,  (34) 
through  (40),  must  be  supplemented  by  boundary  conditions  at  0  ■  0  and  fl  ■  ir, 
and  continuity  conditions  at  0  ■  0^.  By  virtue  of  symmetry,  the  following  con¬ 
ditions  hold  at  0  ■  0. 


U2" 

0 

(41) 

'12  * 

0 

(42) 

ui- 

0 

(43) 

Equations  (30)  and  (31)  together  with  the  conditions  that  and  Oj 2  vani*h 
it  9  *  1  yield  the  following  relations 
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(Bj  -  2)  pB  -  2p  (1  -  B2)1/2  C  -  0  (44) 

2p  (1  -  b£)1/2  A  +  (B2  -  2)p  D  -  0  (45) 

Thus,  the  displacements  and  the  stresses  In  the  elastic  unloading  zone  can  be 

expressed  in  terms  of  tvo  unknown  constants,  say  A  and  B. 

Now,  consider  the  conditions  at  A  m  6  .  An  obvious  condition  follows 

P 

by  comparison  of  (7)  and  (38)  as 


s  -  p  -  1 


(46) 


Since  0  ■  0^  separates  the  plastic  loading  zone  from  the  elastic  unloading  zone, 

e 

a  must  vanish  there.  This  Implies  that  vanishes  at  6  -  6  ;  i.e., 
e  e  p 


-s  I  cos8  +  E'  sln6  ■  0 
e  e 


(47) 


Continuity  of  traction-rates  across  6  -  8^  requires  that 


[Eel  -  0  and  [Eer]  -  0 


(48) 


Here  the  following  notation  has  been  used 


I  ]  -  llm(  )  -  11m (  ) 


e  -►  e  e  -*■  e 

p  p 


Displacement-rates  are  also  continuous  at  6  *  8p.  Consequently, 


(49) 


[Ur]  -  [U0]  -  [Cj]  -  [u2]  =  0 


(50) 


Since  (U,]  •  0,  we  have  [erJ  •  0.  Then,  since  Eg  ■ 


0  and  [Eg]  -  0,  it  follows 
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from  (5)  and  (6)  Chat  [Zr]  *  0.  This  result,  together  with  (48),  Implies 
[d^l  •  0.  From  this  It  follows  that  “  0.  Since  the  attain  rates  and 

the  displacement-rates  are  continuous,  it  follows  from  (11)  through  (13)  that 

lUj]  -  [U£]  -  0  (51) 


at  0  •  8  . 

P 

Since  explicit  expressions  for  the  fields  In  the  elastic  unloading 
zone  are  available,  the  continuity  conditions  (50)  and  (51)  can  be  used  to 
generate  boundary  conditions  for  the  domain  0  <,  0  <_  0^.  The  two  unknown  con¬ 
stants  in  the  unloading  zone,  A  and  B,  can  be  expressed  in  terms  of  U^( 0^  ) 
and  (6p  )  by  the  use  of  (49).  Next,  the  quantities  (U®1)  and  (U^),  which 
are  now  in  terms  of  U. (0  ),  U_(0  ),  s,  6  ,  8.  and  8_  are  computed  and  the 

X  p  L  p  p  L  a 

continuity  conditions  (50)  Imposed.  This  yields  general  relations  of  the  forms 

6i<eP'>  ■  ”i(V  +  *i2(,'V  VV'  (52) 

i2(9p')  '  V#p‘>  +  *22<*-  V  62<  v>  <53) 

The  functions  a  (s,0  ),  which  also  depend  on  8.  and  8T,  are  rather  lengthy,  and 
TO  P  LI 

they  are  not  reproduced  here. 

Numerical  Solution 

Equations  (47),  (52),  and  (53)  provide  a  set  of  boundary  conditions 
for  the  domain  0  <_  8  <_  0^.  The  governing  equations  for  this  domain  are  given 
by  (19)  through  (24).  The  boundary  conditions  at  0  •  0  are  given  by  (41) 
through  (43).  The  unknowns  are  the  functions  depending  on  0,  such  as  U^(0), 
etc,  as  well  as  the  values  of  0^  and  s.  The  problem  defined  In  this  manner  is 
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like  •  nonlinear  eigenvalue  problem,  similar  to,  but  more  complicated  than 
the  one  discussed  for  the  anti-plane  case  by  Achenbach  and  Kanninen  (1978). 

The  problem  must  be  solved  numerically  as  discussed  in  this  section. 

When  (16)  is  used  to  eliminate  in  (19)  through  (24)  and  (47), 
the  latter  can  be  written  in  matrix  form  as 

?(s1;j)  y  +  Jy-0  (54) 

where 

,T  -  [Ur  Dj,  £n,  IJ2.  I12.  I33]  (55) 

in  which  the  superscript  T  denotes  the  transpose.  The  elements  of  the  sym¬ 
metric,  square  matrices  and  K  are  given  in  the  Appendix. 

In  preparation  for  integrating  (54)  numerically,  it  is  convenient  to 
normalize  the  6-variations  by  taking  Ie(0)  »  1  and  defining 

q  -  rn(0)/I22(0)  (56) 

The  nonlinear  two-point  boundary  value  problem  defined  by  (40)  through 
(43),  (50)  and  (54)  through  (56)  is  solved  by  the  shooting  method.  For  pre¬ 
scribed  values  of  a  end  6  and  assumed  values  of  s  and  q,  (54)  is  numerically 
Integrated  over  6  >  0  until  (47)  is  satisfied  at  a  certain  6p.  A  check  is  made 
to  determine  whether  or  not  the  continuity  conditions.  Equation  (51),  are  satis¬ 
fied  for  this  6p.  If  not,  Newton’s  method  is^used  to  establish  new  values  of 
s  and  q  and  the  procedure  is  repeated  until  (51)  is  satisfied. 

A  predictor-corrector  method  was  used  to  integrate  Equation  (54) . 

If  y^  denotes  the  solution  at  0  ■  149,  then  y^+2  at  8  ■  (i+l)66  can  be  written 

?i+i  "  h  +  Zi a®2/2 


(57) 
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and 

?1+1  "  Z'i  +  5  Ae  •  (5®) 

The  coefficient  vector  c  Is  determined  by  numerical  Iterations;  i.e.,  for  an 
assumed  g  (57)  is  substituted  into  (54)  to  determine  y^+^  and  a  new  estimate 
for  c  is  found  from  (58).  This  procedure  is  repeated  until  y^+^  agrees  with 
its  previous  iterant  to  six  significant  figures.  Typically  two  iterations  were 
sufficient  and  the  method  proved  to  be  stable  for  a  as  smell  as  0.005  whereas 
the  finite  difference  technique  used  by  Achenbach  and  Kanninen  (1978)  was  found 
to  be  unstable  for  a  <  0.1.  To  commence  the  numerical  integration  yQ  and  y^ 
are  required.  These  can  be  obtained  by  Introducing  Taylor  series  expansions 
about  6  ■  0  for  Dt  R  and  y  into  (54)  and  equating  coefficients  of  0^  and  0^  to 
xero.  For  the  boundary  conditions  (41)  through  (43)  and  condition  (56),  yQ  end 
y'o  are  given  in  Appendix. 


RESULTS  AND  DISCUSSION 


The  numerical  solution  was  checked  by  comparing  it  with  the  known 
elastic  solution  and  the  quasi-static,  linear  strain-hardening  solution  of 
Asazlgo  and  Hutchinson  (1977).  For  this  comparison  it  is  not  possible  to 
consider  numerically  the  limiting  case  6*0  instead,  6  «  0.001  was  used. 

The  computed  values  of  s  and  0p  were  found  to  agree  with  Amazlgo  and 
Hutchinaon's  values  through  four  significant  figures  for  all  of  the  values 
of  o  that  they  reported;  i.e.  from  a  -  .005  to  1.0. 

For  a  given  value  of  o,  there  exists  a  limiting  crack  speed  above 
which  the  numerical  integration  algorithm  failed  to  converge  to  a  non¬ 
positive  value  of  s.  This  occurred  whenever  the  crack  speed  was  greater 
than  the  Rayleigh  wave  speed  based  upon  the  tangent  modulus;  i.e. 
whenever 

1/2 

3  >  0.57  a  ;  v  -  0.3  .  (59) 

where  0  is  given  by  Equation  (21). 

Numerical  values  of  s  and  6  for  the  condition  of  plane  stress  are 

P 

given  in  Table  1  for  selected  values  of  a  and  0  and  for  v  »  0.3.  Plots  of  s 
versus  a  appear  in  Figure  3.  It  is  apparent  from  this  figure  that  if  the  crack 
speed  is  less  than  approximately  one-half  of  the  Rayleigh  wave  speed  based  upon 
the  tangent  modulus;  i.e.,  if 

1/2 

0  <  0.3  o'  ,  (60) 

the  order  of  the  stress  and  strain  singularity  does  not  differ  significantly 
from  that  for  the  quasi-static  condition.  The  stress  distributions  in  the 
loading  region  for  0  ■  0.25  and  a  -  0.3  are  shown  in  Figure  4.  These  may  be 
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compared  with  the  quasi-static  distributions  for  o  -  0.3  depicted  in  Figure  5. 
While  the  differences  are  not  great,  the  largest  difference  occurs  for  1^. 

Numerical  values  of  s  and  6  are  summarized  in  Table  2  and  plots  of 

P 

a  versus  a  appear  in  Figure  6  for  conditions  of  plane  strain  and  v  •  0.3.  A 
comparison  of  Figures  3  and  6  reveals  that  inertia  effects  have  a  more  pro¬ 
nounced  Influence  upon  the  order  of  the  singularity  for  plane  strain  than  they 
do  for  plane  stress.  The  plane  strain  distribution  of  stresses  in  the  loading 
region  are  shown  in  Figures  7  and  8.  Here  differences  between  the  quasi-static 
distributions  of  Figure  7  and  the  dynamic  counterparts  of  Figure  8  are  more 
perceptible.  In  both  instances  the  high  trlaxlallty  in  the  loading  region  is 
readily  apparent. 
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COHCLPSIONS 

A  basis  for  the  analysis  of  rapid  crack  propagation  taking  direct 
account  of  crack-tip  plasticity  has  bean  provided  In  this  work.  Explicit 
determination  of  the  crack  tip  singularity  and  the  angular  position  at  which 
unloading  takes  place  has  been  sade  for  both  plane  stress  and  plane  strain  con¬ 
ditions.  It  was  found  that,  in  both  conditions,  the  results  are  much  wore  sen¬ 
sitive  to  the  ratio  of  the  slopes  of  the  two  portions  of  the  bilinear  stress- 
strain  relation  used  then  to  the  crack  speed.  This  Indicates  that,  for  a 
singular  finite  element  dynamic  crack  propagation  solution  procedure,  only 
modest  adjustments  In  the  crack  tip  element  formulation  will  be  necessary  to 
accommodate  changes  In  the  crack  speed. 


i 
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APPENDIX 

The  nonzero  elements  of  the  symmetric  metrics  D  (S^)  end  g  for 
plene  strain  ere  es  follows.  Note,  for  plane  stress,  the  lest  column  in  g 
and  0  (Sjj)  end  the  lest  row  in  g,  y,  yQ,  and  D  (S^)  ere  deleted. 

Dn  "  d22  "  ®2,ine  (Al) 


D31- 


-  Dj2  ■  sine 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 


A-  2 


ZABLE  1.  PLANE  STRESS  RESULTS 


o\6 

0.1 

0.25 

0.40 

0.5 

• 

1.0 

-0.500 

-0.500 

-0.500 

-0.500 

0,7 

-0.460 

-0.457 

-0.439 

0.5 

-0.419 

-0.411 

-0.371 

0.3 

-0.355 

-0.339 

0.2 

-0.306 

-0.282 

0.1 

-0.232 

e 

P 

1.0 

1.408 

1.496 

1.706 

0.750 

0.7 

1.426 

1.516 

1.712 

0.5 

1.425 

1.517 

1.719 

0.3 

1.401 

1.498 

0.2 

1.370 

1.475 

0.1 

1.304 

1.0 

-0.500 

-0.500 

-0.500 

0.7 

-0.472 

-0.468 

-0.445 

0.5 

-0.440 

-0.426 

-0.336 

0.3 

-0.369 

-0.320 

0.2 

0.1 

-0.292 

-0.180 

-0.161 

1.0 

1.550 

1.630 

1.771 

0.7 

1.624 

1.699 

1.798 

0.5 

1.718 

1.768 

1.834 

0.3 

1.871 

1.847 

0.2 

1.967 

1.802 

0.1  2.091 


FIGURE  1.  CRACK  TIP  GEOMETRY 


FIGURE  2.  BILINEAR  STRESS-STRAIN  CURVE 


FIGURE  3.  ORDER  OF  THE  CRACK  TIP  SINGULARITY  IN  PLANE  STRESS 
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FIGURE  5.  QUASI-STATIC  CRACK  TIP  STRESSES  IN  PLANE  STRESS 


FIGURE  6.  ORDER  OF  THE  CRACK  TIP  SINGULARITY  IN  PLANE  STRAIN 
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FIGURE  7.  DYNAMIC  CRACK  TIP  STRESSES  IN  PLANE  STRAIN 


